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We report MC calculations of perturbative coefficients for lattice scalar field theory in dimensions 1, 2 and 
3, where the large field contributions are cutoff. This produces converging (instead of asymptotic) perturbative 
series. We discuss the statistical errors and the lattice effects and show that accurate calculations are possible 
even in a crossover region where no approximation works. We show that the field cutoff is also a UV regulator. 
We point out the relevance for QCD questions discussed by Tomboulis and Trottier at this conference. 



A central problem in quantum field theory is 
the presence of different behaviors at different 
scales. This is clearly the case for QCD pQ 
where the short distance behavior can be de- 
scribed using perturbation theory, while confine- 
ment and other non-perturbative phenomena ap- 
pear at large distance. A somehow similar situ- 
ation is encountered in scalar and spin models. 
We limit here the discussion to the symmetric 
phase where a trivial high-temperature (strong 
coupling) fixed point is present. Perturation the- 
ory generically fails to provide a proper descrip- 
tion of the RG flows near the HT fixed point. 
From this point of view, the fact that the pertur- 
bative series are asymptotic is not surprising. 

Asymptotic series can be a serious problem for 
problems where the strong interactions correction 
are taken into account. A good example is the 
hadronic width of the Z , where the third order 
corrections are about two-third of the second 
order corrections and larger than the experimen- 
tal error bar for the combined LEP experiments. 
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Perturbative methods also play an important role 
in the improvement methods that were used 
to obtain unprecedented accuracy for quantities 
that can be compared with experiments @]. Given 
these recent successes, the standards of accuracy 
are raised and it is likely that the lack of conver- 
gence of QCD series will need to be addressed in 
the near future. 

A simple wav |5!6) to convert an asymptotic se- 
ries into a convergent one consists in introduc- 
ing a large field cutoff. It has been shown with 
three (ft 4 problems t 6 that the modified series con- 
verge toward values which are exponentially close 
to the exact ones. At a fixed order, it is possi- 
ble to chose the field cutoff in order to optimize 
the accuracy [7|. The three examples discussed in 
Ref. :6. can be solved accurately with numeri- 
cal methods |8l9j . However, for generic problems, 
this kind of calculation can be quite difficult, es- 
pecially in a crossover region (illustrated in Fig. 
3) where neither semi-classical methods or HT ex- 
pansion are available. This difficulty reflects the 
fact that, in general, the interpolation between 
RG fixed points is a difficult non-linear problem 
that has only be solved for simplified models (see 
[TTJ] for an example). 

When no other methods are available, one has 
to resort to the MC method to perform calcula- 
tions of the modified perturbative coefficients. In 
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the following, we consider lattice scalar models 
with nearest neighbor interactions and one com- 
ponent at each site, in 1, 2 or 3 dimensions. For a 
lattice with N s sites, we will consider the follow- 
ing quantities relevant for perturbation theory: 
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M = (<(£^) 2 >-<£<^> 2 )/^- (3) 

X X 

The average <> correspond to a Gaussian model 
where the integration at each site goes from 
—4>max to 4>max- In D = 1, the continuum val- 
ues for 4> m ax = oo are Aq = 1/2, Ai = 3/4 and 
A 2 = 21/8. For D=2 or 3, A is UV divergent. 
Two kinds of errors that should be taken into ac- 
count: the statistical errors and the errors due to 
the finite lattice spacing. 



we can use Ref. ^2] to estimate the statistical 
error for correlated data. 
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Figure 2. Illustration of the power corrections to 
A\ in D = 1. We have subtracted the estimated 
continuum value c. 
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Figure 1. Correlation functions for the observable 
Aq in D — 3 as a function of the separation n in 
the configuration number. 



Second, we need to take into account the lat- 
tice effects. For small enough lattice spacing a, a 
simple power behavior c + a a is observed, c being 
the continuum value that we are trying to calcu- 
late. The power behavior can be extracted from 
the statistical fluctuations (which increase when 
a becomes small!) for an intermediate range of 
lattice spacing, using a nonlinear fit. This is il- 
lustrated in Fig. 2. More detail will be provided 

inna. 

The comparison between the MC results ob- 
tained with the method described above and the 
accurate numerical results obtained with meth- 
ods inspired by Ref. [5] is shown below in the 
case D = 1 . The MC estimates for A\ for var- 
ious 4> m ax are compared with accurate values in 
the table below 



The statistical errors decrease as the inverse 
square root of the number of decorrelated con- 
figurations. The correlation times are estimated 
from the exponential decay of the correlations. If 
we denote by Oi the value of the observable O in 
a configuration i, the correlations C(n) are sub- 
tracted averages of OiOi+ n . Fig. 1 illustrates this 
exponential decay in an example. Alternatively, 



(pin ax 


MC 


Num. 


3 


0.738(4) 


0.7333 


2.5 


0.646(12) 


0.6482 


2 


0.422(9) 


0.4315 


1.5 


0.178(2) 


0.1807 


1 


0.03987(3) 


0.03995 



The results for A\ and A 2 are shown in Fig. 3. 
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The MC estimates for A 2 require a subtraction 
and are less accurate. We expect that better re- 
sults can be obtained by improving the statistics 
and refining the nonlinear fit method. 
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Figure 3. MC (circles) versus accurate numeri- 
cal results (continuous line) for A\ (left) and A 2 
(right). All the values have been divided by their 
respective infinite cut limits 3/4 and 21/8. 



It is possible to calculate approximately the 
values at the bottom (top) of the curves of Fig. 
3 using a HT expansion (semi-classical methods) 
and reliable interpolation methods are being de- 
veloped to apprach the crossover region. 

Aq diverges like Ln(a) in D = 2 and 1/a in 
D = 3, when the lattice spacing a — > 0. However, 
the field cut takes care of this divergence. This 
is illustrated in Fig. 4 for D — 2. This regu- 
larization has a simple explanation 13 : since the 
average (1) provides a positive measure, we can 
obtain a bound by replacing <p^. by its maximal 
value. This yields the bound A < ^ax- The 
field cut can be thought as produced by an inter- 
action of the form (<j) / 4> m ax) s in the large s limit. 
Consequently, this should not affect the universal 
features of the model and it should be possible 
to calculate the critical exponents using modified 
pcrturbative methods. We are planning to extend 
these methods to gauge theories. 
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Figure 4. Aq for D = 2 as a function of the lattice 
spacing for various field cuts. The continuous line 
is the (calculable) lattice result (without cut). 
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